Mon. Not. R. Astron. Soc. 000, [TV?? (2012) Printed 30 July 2012 (MN ET E X style file v2.2) 



Pulse-driven nonlinear Alfven waves and their role in the 
spectral line broadening 

P. Chmielewski, 1 A.K. Srivastava, 2 K. Murawski, 1 and Z.E. Musielak 3 ' 4 

1 Group of Astrophysics, UMCS, ul. Radziszewskiego 10, 20-031 Lublin, Poland 

2 Aryabhatta Research Institute of Observational Sciences (ARIES), Manora Peak, Nainital-263 129, Uttarakhand, India 

3 Department of Physics, University of Texas at Arlington, Arlington, TX 76019, USA 

4 Kiepenheuer-Institut fur Sonnenphysik, Schdneckstr. 6, 79104 Freiburg, Germany 



Released 2002 Xxxxx XX 



ABSTRACT 

We study the impulsively generated non-linear Alfven waves in the solar atmosphere, 
and describe their most likely role in the observed non-thermal broadening of some 
spectral lines in solar coronal holes. We solve numerically the time-dependent mag- 
netohydrodynamic equations to find temporal signatures of large-amplitude Alfven 
waves in the model atmosphere of open and expanding magnetic field configuration, 
with a realistic temperature distribution. We calculate the temporally and spatially 
averaged, instantaneous transversal velocity of non-linear Alfven waves at different 
heights of the model atmosphere, and estimate its contribution to the unresolved non- 
thermal motions caused by the waves. We find that the pulse-driven nonlinear Alfven 
waves with the amplitude A v —50 km s -1 are the most likely candidates for the non- 
thermal broadening of Si VIII A1445.75 A line profiles in the polar coronal hole as 
reported by Banerjee et al. (1998). We also demonstrate that the Alfven waves driven 
by comparatively smaller velocity pulse with its amplitude A v =25 km s^ 1 may con- 
tribute to the spectral line width of the same line at various heights in coronal hole 
without any significant broadening. The main conclusion of this paper is that non- 
linear Alfven waves excited impulsively in the lower solar atmosphere are responsible 
for the observed spectral line broadening in polar coronal holes. This is an important 
result as it allows us to conclude that such large amplitude and pulse-driven Alfven 
waves do indeed exist in solar coronal holes. The existence of these waves and their 
undamped growth may impart the required momentum to accelerate the solar wind. 

Key words: Sun: atmosphere - (magnetohydrodynamics) MHD - waves. 



1 INTRODUCTION 

Alfven waves are difficult to observe and yet several indica- 
tions of their presence in the solar atmosphere were found 
during SOHO and TRACE space missions. More direct evi- 
dence for the existence of these waves in different regions of 
the solar atmosphere was given by high-resolution observa- 
tions carried out by the Solar Optical Telescope (SOT) and 
the X-Ray Telescope (XRT) onboard the Hinode Solar Ob- 
servatory. According to Okamoto & De Pontieu (2011), De 
Pontieu et al. (2007), and Cirtain et al. (2007), the signa- 
ture of Alfven waves were observed respectively in spicules 
and X-ray jets using SOT and XRT instruments. Interpreta- 
tions of these observational results were discussed by Erdelyi 
& Fedun (2007), Tomczyk et al. (2007), and Antolin et al. 
(2009). 

Observational evidence for the existence of torsional 
Alfven waves in the solar atmosphere was also reported by 
Jess et al. (2009), who analyzed Ha observations obtained 



with high spatial resolution by the Swedish Solar Telescope 
(SST). They interpreted the data in terms of Alfven waves 
in the solar chromosphere, with periods from 12 min down 
to the sampling limit of the observations near 2 min, with 
maximum power near 6-7 min. The authors concluded that 
the amount of energy carried by such transversal waves was 
sufficient to heat the solar corona (Dwivedi & Srivastava 
2006). 

These recent discoveries of Alfven waves in the solar at- 
mosphere well-justified extensive studies of these waves that 
have been performed by numerous investigators in the last 
four decades (e.g., Hollweg 1985; Roberts 1991; Musielak & 
Moore 1995; Narain & Ulmschneider 1996; Roberts 2004; 
Antolin & Shibata 2010, and references cited there). The 
studies have covered both linear (e.g., Hollweg & Isenberg 
2007; Murawski & Musielak 2010) and non-linear (e.g., Ver- 
dini & Velli 2007; Verdini et al. 2009; Matsumoto & Shibata 
2010) Alfven waves, and different aspects of their genera- 
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tion, propagation and dissipation have been investigated. 
The specific objectives of these studies were to understand 
the role of Alfven waves in the atmospheric heating and in 
the acceleration of supersonic solar wind. The fast compo- 
nent of the solar wind originates in solar polar coronal holes 
(e.g., Hassler et al. 1999; Tu et al. 2005, and references cited 
there), which are the regions where the non-thermal broad- 
ening of spectral lines has also been observed (e.g., Hassler 
et al. 1990; Banerjee et al. 1998; Moran 2003; O'Shea et al., 
2005; Dolla & Solomon 2008). Similar observations have also 
been done in the equatorial corona (Harrison et al., 2002). 
These authors proposed that the radially propagating Alfven 
waves may result in the non-thermal broadening of spectral 
line widths. 

These observations were also investigated analytically 
by Pekiinlii et al. (2002) and Dwivedi & Srivastava (2006). 
Recently, the observations from Bemporad & Abbo have 
given the first signature of the dissipation of Alfven waves 
neat the solar limb, which confirm the theoretical models of 
Dwivedi & Srivastava (2006), and Srivastava et al. (2007). 
Moreover, Zaqarashvili et al. (2006) have reported that the 
resonant energy conversion from Alfven to acoustic waves 
in the region where plasma /3 approaches unity in the solar 
atmosphere. This conversion can be responsible for the spec- 
tral line width variation. However, this theory only explains 
the most probable cause of the line-width reduction, which 
was observed only by O'Shea et al. (2005) in solar coronal 
hole. Additional problem is the fact that there is not enough 
observational evidence for the resonant energy conversion in 
the solar corona (e.g., Srivastava & Dwivedi 2010; McAteer 
et al. 2003). Hence, new studies of the role played by Alfven 
waves in the observed spectral line broadening were neces- 
sary. 

In this paper, we numerically study the behavior of 
large-amplitude (non-linear) Alfven waves in a model that 
resembles a solar coronal hole. Our main objective is to de- 
termine the role played by these waves in the spectral line 
broadening as observed in the coronal holes (e.g., Banerjee 
et al. 1998; Dolla & Solomon 2008). We find an agreement 
between our numerical results of non-linear Alfven waves 
and the computed line broadening of constituted synthetic 
spectra at different heights in model coronal hole and the 
observational data. This allows us to conclude that large- 
amplitude Alfven waves are responsible for the observed 
non-thermal broadening of the spectral lines in the coro- 
nal holes. Our result is important because it is an indirect 
evidence for the existence of non-linear Alfven waves in solar 
coronal holes. 

The outline of the paper is as follows: our numerical 
model is described in Sec. 2; a brief description of the used 
numerical code and the form of initial perturbations are 
given in Sec. 3; the results of our numerical simulations are 
presented in Sec. 4; comparison of our results to the ob- 
servational data is given in Sec. 5; the obtained results are 
discussed in Sec. 6; and our conclusions are given in Sec. 7. 



2 NUMERICAL MODEL OF ALFVEN WAVES 
2.1 MHD equations 

Our model of the solar atmosphere contains a 
gravitationally-stratified magneto-plasma, which is de- 



scribed by the following set of ideal magnetohydrodynamic 
(MHD) equations: 

do 



dt 



+ V ■ (qV) = , 



(1) 



e ^X + e CV.V)V=-Vp+-(VxB)xB + gg, (2) 



dB 

~dt 



dt 



V x (V x B) , (3) 
V • B = , (4) 
+ V ■ Vp = - 7 pV ■ V , (5) 



P = — QT. 

m 



(6) 



Here q is mass density, V and B are vectors of respectively 
the flow velocity and the magnetic field, p is gas pressure, 
7 = 5/3 is the adiabatic index, g = (0, — <7,0) is a vector 
of gravitational acceleration with its value g = 274 m s -2 , 
T is a temperature, m is a mean particle mass and fee is a 
Boltzmann's constant. 



2.2 A model of the solar atmosphere 

We consider a model of the solar atmosphere with an invari- 
ant coordinate (d/dz = 0) and allow the z-components of 
velocity (V z ) and magnetic field (B z ) to vary with x and y. 
In such 2.5D model, the solar atmosphere is in static equilib- 
rium (V e = 0) with force- and current-free magnetic field, 
i.e., 



(VxB e )xB c = 0, V x B c 



0. 



(7) 



Henceforth the subscript e corresponds to equilibrium quan- 
tities. 

In our model of the atmosphere, a curved magnetic field 
is given by 



B, 



VxA c 



where the magnetic flux function (A c ) has the form 



Ab-Bq cos 



(aI) 6XP (■ 



A B 



(8) 



(9) 



Here z is a unit vector along the z-direction and Bo is the 
magnetic field at the reference level, y = y r , that is chosen 
at y t — 10 Mm. We set and hold fixed Bp in such a way 
that the Alfven speed, ca = Bo/ \J /j,g e (y = y T ) is ten times 

higher than the sound speed, c s = yj jp c {y = J/r) / Q c (y = y r ). 
Such a choice of B e results in Eq. (O being satisfied. Here 
As = 2L/tt denotes the magnetic scale-height, and L is 
a half of the magnetic arcade width. As we aim to model 
a polar region, we take L = 75 Mm and keep it fixed in 
our calculations. For such a choice, the magnetic field lines 
are weakly curved (not shown) and represent the open and 
expanding field lines similar to the coronal holes. 

As a result of Eq. (J7J), the pressure gradient is balanced 
by the gravity force 

-Vp c + (? c g = 0. (10) 

Using the ideal gas law of Eq. (|6]) and the y-component of 
the hydrostatic pressure balance indicated by Eq. (|10|l . we 
express equilibrium gas pressure and mass density as 



Pe(y) = Po exp 



dy 

w 



Qe(y) = 



Pc(y) 
<?A(y) 



(11) 
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Figure 1. Equilibrium profile of the temperature (top) and the 
Alfven speed (bottom). 



Here 



A(») = 



fc B r c ( y ) 

mg 



(12) 



is the pressure scale-height, and po denotes the gas pressure 
at the reference level. 

We adopt a realistic model of the plasma temperature 
profile (Vernazza et al. 1976), displayed in Fig.[T], top panel. 
Temperature attains a value of about 6 x 10 3 K at y = 1.5 
Mm and it grows up to about 1.5 x 10 6 K in the solar corona 
at y = 10 Mm. Higher up temperature is assumed constant; 
note that this profile is more realistic than the one used 
by Murawski & Musielak (2010). The temperature profile 
determines uniquely the equilibrium mass density and gas 
pressure profiles. Both g e (y) and p e (y) experience a sudden 
drop at the transition region that is located at y ~ 2.7 Mm 
(not shown). 

In this model the Alfven speed, ca, varies only with y 
and is expressed as follows: 



B e~ 



(13) 



Its profile is displayed in Fig. [TJ bottom panel. Note that 
the Alfven speed in the chromosphere, c\(y = 1.75 Mm), is 
about 25 km s~ . The Alfven speed rises abruptly through 



the transition region reaching a value of cx{y = 10 Mm) = 
10 3 km s" 1 (Fig. [TJ bottom). The growth of c A (y) with 
height results from a faster decrement of Q c {y) than B c (y) 
with the height. 

The realistic solar atmosphere above the polar corona 
reveals complexity of its plasma and magnetic field struc- 
ture. The magnetic field configuration in the polar coronal 
hole can be approximated by expanding coronal funnels in 
the lower part of the atmosphere and comparatively smooth 
open field lines in its upper part (Banerjee et al. 1998; Hack- 
enberg et al. 2000). This implies that the field structure and 
magnetic scale-height vary above the polar corona. More- 
over, the field configuration changes from dipolar to multi- 
polar during the transition from the solar minimum to its 
maximum. Despite these well-known variations, we assume 
that the magnetic field scale-height is fixed at a resonable 
value of L = 75 Mm in our simulation domain, which rep- 
resents the weakly curved and open field lines of the polar 
coronal hole for the quiet minimum phase of the Sun. We 
want to point out that this assumption does not affect the 
validity of our numerical results. 



3 NUMERICAL SOLUTIONS OF MHD 
EQUATIONS 

Equations (HJ)-© are solved numerically with a use of the 
FLASH code (Fryxell et al. 2000; Lee & Deane 2009; Lee 
et al. 2009). This code implements a second-order unsplit 
Godunov solver with various slope limiters and Riemann 
solvers as well as Adaptive Mesh Refinement (AMR) (Mac- 
Neice et al. 1999). We use the minmod slope limiter and the 
Roe Riemann solver (e.g., Toro 2009). We set the simulation 
box as (—5 Mm, 5 Mm) x ( — 1 Mm, 84 Mm) and impose fixed 
in time boundary conditions for all plasma quantities in the 
x- and y-directions, while all plasma quantities remain in- 
variant along the z-direction. In our present work, we use 
AMR grid with a minimum (maximum) level of refinement 
set to 3 (8). The refinement strategy is based on control- 
ling the numerical errors in mass density. Blocks are denser 
below 3 Mm and vertically along the region of Alfven wave 
propagation. Every numerical block consists of 8 x 8 iden- 
tical numerical cells. This results in an excellent resolution 
of vital spatial profiles and greatly reduces the numerical 
diffusion at these locations. 



3.1 Initial perturbations 

We perturb initially (at t — s) the model equilibrium, de- 
scribed in Sec. 12.21 by a Gaussian pulse in the z-component 
of velocity given by 



V z (x, y, t — 0) = A v exp 



(a; - x ) 2 + (y- y ) 2 



(14) 



where A v is the amplitude of the pulse, (xo,yo) is its ini- 
tial position and w denotes its width. We set w = 1 Mm, 
(xo = 0, yo = 1.75) Mm and consider two cases: (a) A v = 25 
km s _1 ; (b) A v = 50 km s _1 . 

Note that in the 2.5D model, we developed, the Alfven 
wave decouples from magnetoacoustic waves and it can be 
described by V z (x, y, t). As a result, the initial pulse triggers 
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0.005 



0.000 



for A v — 50 km s -1 collected at (x = 0, y = 49) Mm are 
shown in right-bottom panel of Fig. [3] Since the wave is non- 
linear, it generates a vertical flow (Fig. [3] middle-bottom 
panel) and mass density perturbations (Fig. [3] left-bottom 
panel), which are driven by the ponderomotive force. 

We collect the temporally and spatially averaged 
transversal velocity component V z at each height of the sim- 
ulation domain (y ma m = 84 Mm) generated by the pulses of 
different wave amplitudes, A v — 25 and 50 km s _1 (see Fig. 
3). If we assume that the y-direction of the simulation box 
is placed along the outward open magnetic field of a polar 
coronal hole, then at each height the transversal amplitude 
of Alfvenic perturbations may contribute to the otherwise 
unresolved non-thermal motions by affecting the line-width 
of the observed spectral lines. Finite-amplitude Alfven waves 
propagating along the open magnetic field lines of the polar 
coronal holes can perturb the plasma velocity, which causes 
positive and negative Doppler shifts that can be detected 
as a line width broadening or line width variation (Baner- 
jee et al. 1998; Harrison et al. 2002; Moran 2003; O'Shea et 
al. 2005; Dwivedi & Srivastava 2008). We now compare our 
results of numerical simulations with the available observa- 
tional data. 



Figure 2. Transverse velocity V z profiles at t = 20 s and t = 120 
s for A v = 25 km s . 



Alfven waves that in the linear limit are described by the 
wave equation 



9 2 K 2 , x d 2 V* 



dt 2 



dy 2 



(15) 



4 RESULTS OF NUMERICAL SIMULATIONS 

We simulate impulsively excited non-linear Alfven waves 
and investigate their propagation along the open magnetic 
field lines of a coronal hole model in the outward direction. 
It should be noted that the effect of the inhomogeneities 
across the magnetic field lines are not included in our ap- 
proach. The waves are generated by the transversal velocity 
pulse perpendicular to the magnetic isosurface (X-Y) in the 
z-direction. This pulse is described by Eq. 1)140 , 

First, we consider the pulse amplitude equal to 25 
km s _1 . Spatial profiles of V z for A v =25 km s _1 and the cor- 
responding time signatures obtained in our numerical simu- 
lations of Alfven waves are shown respectively in Figs. [2] and 
[21 top-right panel. Perturbations of K propagate essentially 
along magnetic field lines, which is displayed on spatial pro- 
files of V z at t = 20 s and 120 s (Fig. [2j . A part of the whole 
simulation region is illustrated. Note that the Alfven wave 
experiences an acceleration at that height as CA(y) increases 
up to 1 Mm s _1 at the transition region (Fig.[T] right). This 
acceleration of Alfven wave is caused by significant fall in 
mass density at essentially constant magnetic field. Above, 
the Alfven wave propagates with slightly increased ampli- 
tude and almost constant velocity. The spatial resolution of 
the transition region is set to be Ax = Ay = 0.04 Mm in 
the simulation domain. 

We also simulate the Alfven wave using the large 
transversal pulse of A v — 50 km s _1 that generates the large- 
amplitude, non-linear Alfven waves. Time-signatures of V z 



5 COMPARISON TO OBSERVATIONAL DATA 
5.1 Line width broadening 

The non-thermal motion is a prominent candidate that may 
modify the observed line profiles by its contribution to the 
full width at half maximum (FWHM) of the spectral lines. 
Banerjee et al. (1998) have estimated the radial variation 
of the non-thermal velocity in the polar coronal hole by de- 
ducing its contribution to the observed line profiles of Si 
VIII A1445 A . The sophisticated line-width measurements 
and thus the observed variation of non-thermal velocity are 
very few in the available literature. Using SOHO/CDS ob- 
servations, O'Shea et al. (2005) have shown the line width 
increment in the inner corona and then its decrement be- 
yond radial distance 1.21 R in coronal hole. Here R is the 
radius of the Sun. However, due to the unavailability of the 
instrumental width of SOHO/CDS slits, the exact informa- 
tions about the non-thermal velocity could not be derived, 
so the estimate gave only the uncorrected line-width varia- 
tion with the height. Moreover, Banerjee et al. (1998) have 
approximately re-produced the trend of the spatial variation 
of non-thermal velocity in the polar coronal hole similar to 
Banerjee et al. (1998) using 2" slit of Hinode/EIS. Note that 
the exact instrumental width of EIS slits were variables with 
wavelength of the observed spectral lines. 

By using the temperature averaged over the simulation 
domain, the averaged density generated from the simula- 
tions, and the temporally and spatially averaged transversal 
speed (V z ) of Alfven wave generated by the velocity pulse A v 
= 50 km s - , we have estimated the equivalent FWHM (a) 
of its 1445.75 A line at different heights using the relation 
(Mariska 1992) 



41n2(^ 2 f^ + f 



i 2 



(16) 



It should be noted that the averaged wave velocity am- 
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Figure 3. Results of the numerical simulation of pulse-excited Alfven waves: time-signatures of p, V y and V z for A v = 25 km s _1 (top 
panels) and A v = 50 km s _1 (bottom panels) collected at (x = 0, y = 49) Mm. 



plitude can be scaled appropriately in terms of non-thermal 
speed as £ 2 = 0.5 V 2 z by taking into consideration of polar- 
ization and direction of the propagation of wave w.r.t line- 
of-sight (Banerjee et al. 1998). We adopt the same scaling 
throughout the paper, although different authors use dif- 
ferent scalings based on their assumption of the degree of 
freedom of wave motion (Doha & Solomon 2008). At any 
particular height in the model coronal hole, we average the 
wave velocity amplitude of transversal oscillation in tempo- 
ral domain between i a -250 to £ a +250 s, while in the spatial 
domain over the entire pulse width, where £ a is the arrival 
time of a wave signal to the detection point. Therefore, the 
resultant transversal motion contributes to the non-thermal 
unresolved motion at each height in the model coronal hole 
plasma as observed in form of spectral line broadening by 
various authors, e.g., Banerjee et al. (1998) and Dolla & 
Solomon (2008). 

5.2 Line profiles obtained with A v = 50 km s _1 

We computed the synthetic line profiles (cf., Fig. 4) of Si 
VIII A1445.75 A at three different heights y = 20, 40, and 70 
Mm based on the inclusion of the non-thermal contribution 
of Alfven wave excited by a pulse with its amplitude A v — 50 
km s - . The CHIANTI atomic data base is used to produce 
the synthetic line profiles by computing the averaged density 
(p), temperature (T), and FWHM estimated by averaged ve- 
locity (V*). The CHIANTI SSW routine "synthetic.pro" is 
used to produce the line profiles; note that the routine takes 
into account the ionic equilibrium as reported by Mazzotta 
et al. (1998), the coronal hole differential emission measure 
(DEM) values, and coronal abundances as available in CHI- 
ANTI atomic data base. Although, the major input param- 
eters are the averaged density (p) (or pressure), tempera- 
ture (T), and FWHM (see Eq. [T6)) derived from the simula- 



tion, while other parameters mentioned above are considered 
as optional input parameters in the calculation of synthetic 
spectrum. 

In the present case, we compare the deduced theoreti- 
cal line width with the corrected observed one as reported 
by Banerjee et al. (1998) using SoHO/SUMER. There- 
fore, we neglected the implication of instrumental width of 
SoHO/SUMER slits in Eq. (16). We synthesized the line pro- 
files of the Si VIII 1445.75 A line but did not calculate the 
continuum as we wanted to only show the shape of the line 
profiles at various heights and their broadening that com- 
pares well with the observations of Banerjee et al. (1998). 
The line width of synthetic line profiles of Si VIII 1445.75 
A at 20 Mm, 40 Mm, and 70 Mm (Eq. 16) are respectively 
353 mA , 358 mA , and 363 mA when we consider the non- 
thermal contribution of impulsively excited Alfven wave by 
a pulse with its amplitude ^4 V =50 km s _1 . It is clear from 
Table 2 of Banerjee et al. (1998) that the estimated cor- 
rected line-widths of Si VIII A1445.75 A at heights 27 arc- 
sees (~19.6 Mm), 57 arcsecs (~ 41.3 Mm), 98 arcsecs (~71 
Mm) are respectively ~292 mA , ~334A , and ~369 mA . 

We also adopted and tested a different way of calculat- 
ing the line-profiles at a particular height (e.g., 15 Mm) in 
the model coronal hole. Using the results of our numerical 
simulation, we derived the plasma parameters (p, T) as well 
as the transversal speed V z for each instant set (X, Y=fixed, 
t) in which X-space is averaged over 1.0 Mm spatial scale. 
The time instant was considered with each 10 s resolution 
between £-250 s to £+250 s. We collected the parameters at 
each instant, and then derived the averaged FWHM, which 
was found to be 284 mA for a pulse ^4 V =50 km s _1 . The 
result closely matches the broadened line width of 291 mA 
as reported by Banerjee et al. (1998) at a height of 23 arc- 
sec. Therefore, the alternative estimations are also consistent 
with our findings. 
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Although our theoretically estimated line-width also 
shows increment with height in our model coronal hole, the 
spatial gradient of increment is rather flat compared to the 
observations by Banerjee et al. (f998). The most likely rea- 
son is that we considered the impulsive excitation of the 
pulse-driven Alfven wave that also contributes more to the 
un-resolved non-thermal motion near the wave source region 
in the lower solar atmosphere. Therefore, the computed line- 
width, ,for instance at 27 arcsecs (—19.6 Mm), is higher than 
the observed line width given by Banerjee et al. (1998) at 
the same height. However, the theoretically estimated and 
the observed line- widths at coronal heights e.g., at 40 and 
70 Mm, closely match each other. 

The results presented in Fig. 4 clearly show that the 
synthetic line profiles computed with the wave amplitude 
A v — 50 km s _1 closely approximately mimic the observed 
spectrum of Si VIII A1445.75 A obtained by the detector-B 
on SOHO/SUMER and reported by Banerjee et al. (1998). 
This provides anevidence that the large-amplitude, pulse- 
excited Alfven wave can contribute upto some extent into 
the broadening of the spectral lines observed by Banerjee et 
al. (1998). 

5.3 Line profiles obtained with A v = 25 km s _1 

We again derived the synthetic line profiles (cf., Fig. 5) of Si 
VIII A1445.75 A (Doha & Solomon 2008) at three different 
heights y = 43.5, 65.25, and 72.5 Mm, based on the inclusion 
of the non-thermal contribution of Alfven wave excited by 
the pulse with its amplitude A v = 25 km s _1 . The CHIANTI 
atomic data base and the method described above was again 
used to produce the synthetic line profiles. The obtained 
results show that there is not very significant broadening 
of line profiles in our model coronal hole for these smaller 
amplitude Alfven waves; for the observational data see Doha 
& Solomon (2008). 

The line widths of Si VIII A1445.75 A synthetic line pro- 
files at 43.5 Mm (~60"), 65.25 Mm (~90"), and 72.5 Mm 
(—100") (Fig. 5) above the limb in our coronal hole model 
are estimated respectively as —289 mA —290 mA and 291 
mA when we consider the non-thermal contribution of im- 
pulsively excited Alfven wave by a pulse strength of ^4 V =25 
km s _1 . Doha & Solomon (2008) have reported the Gaus- 
sian half line- widths of various coronal spectral lines (e.g., 
Si VIII A1445.75 A Fe XII A1242, Mg X A624 A ) between 
-42 Mm (57") and -100 Mm (140") in the polar corona 
hole (cf., their Fig. 8). This shows that our theoretical esti- 
mations are below the observationally estimated line-width 
as reported by Dolla & Solomon (2008) for Si VIII A1445.75 
A between 40-75 Mm in coronal hole. For example, the es- 
timated observed half line- widths at 1/ yjje) as reported by 
Dolla & Solomon (2008) (cf., their Fig 8), after conversion in 
FWHM of Si VIII A1445.75 A at heights -43.5 Mm (-60"), 
- 65.25 Mm (-90"), and -72.5 Mm (~100"), respectively 
give -362 mA -368A , and -376 mA . 

Our theoretically estimated line- width shows very small 
increment w.r.t. height in the coronal hole model between 
40-75 Mm height. Now, the observed width of Si VIII 
1445.75 A also varies very little but with comparatively 
higher gradient than theoretical values. The most proba- 
ble reason is that the plasma and magnetic field properties 
of the observed coronal hole differs from the assumed con- 



ditions in our model, in which the Alfven wave of smaller 
amplitude propagates without much growth as well as its 
contribution in the line broadening. 

Our numerical results obtained for non-linear Alfven 
waves with the amplitude A v =25 km s _1 show that the re- 
sulting contribution of these waves to non-thermal motions 
approximately qualitatively match the observational data of 
spectral line width as observed by Dolla & Solomon (2008) 
for Si VIII A1445 A line of SoHO/SUMER, as both show the 
lesser spectral line broadening. However, the theoretically 
estimated widths are comparatively smaller than the ob- 
served ones at various heights in coronal hole. Although, the 
theoretically derived line- widths only approximately match 
the observed ones at various heights in the coronal hole, 
there is no clear evidence of significant line broadening in 
these cases. 



6 DISCUSSION 

As already mentioned in Sec. 1, there are several analytical 
models that attempted to explain the observed non-thermal 
broadening of spectral lines in the solar corona by using lin- 
ear Alfven waves (e.g., Harrison et al. 2002 ; Pekiinlii et al. 
2002; Dwivedi & Srivastava 2006). Under the linear theory 
of propagating Alfven waves from the solar photosphere up- 
ward into the solar corona along the magnetic field lines, it 
is found that the non-thermal velocity is inversely propor- 
tional to the quadratic root of the observationally estimated 
electron density that matches well the characteristics of such 
waves (e.g., Banerjee et al. 1998; Moran et al. 2003; Baner- 
jee et al. 2009). The falls-off of the density with height can 
amplify the propagating Alfven wave and increase its group 
velocity that can contribute in the rms wave amplitude and 
thus in the non-thermal spectral line- width broadening with 
height in the solar atmosphere. 

The most relevant numerical work has been done by 
Kudoh & Shibata (1999) , who used numerical simulations 
to investigate the role of torsional Alfven waves in both the 
formation of spicules as well as the observed non-thermal 
broadening of spectral lines in the solar corona. They con- 
sidered Alfven waves to be generated by perturbations of 
1 km s _1 in the solar photosphere. According to these au- 
thors, the waves can lift the spicules up to 5000 km, provide 
the energy flux 3.0 x 10 5 ergs s _1 cm~ 2 , and produce the 
non-thermal broadening of emission lines in the solar corona 
with approximately 20 km s _1 . 

It is clear from the above discussion that neither previ- 
ous analytical methods nor previous numerical simulations 
have been able to explain the full profile of the observed 
spectral line-width variation with height in the polar coro- 
nal holes. Nevertheless, it is noteworthy here that previous 
observational results have been matched with linear propa- 
gating waves coming from photosphere (Banerjee et al. 1998, 
2009). They have showed, in linear regime, that the rms wave 
amplitude is related with some powers of mass density and 
magnetic field variations. In this paper, we are not attempt- 
ing this indirect analytical model explanations but instead 
we directly compute the temporally and spatially averaged 
transversal velocity component, V z , from our pulse-excited 
Alfven wave model and determine its contribution to rms 
wave velocity amplitude. It should be noted then that the 
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Figure 4. The synthetic line profiles of Si VIII 1445.75 A line at three different heights y = 20, 40, and 70 Mm to compare with the 
observations of Banerjcc ct al. (1998), based on the inclusion of the non-thermal contribution of Alfven wave excited by pulse amplitude 
A v = 50 km s — 1 . The line-profile is deduced by considering the ionic equilibrium as reported by Mazzotta et al. (1998), the coronal hole 
DEM values, and coronal abundances as available in CHIANTI atomic data base. The simulated temperature and density are also used 
as input parameters. 



non-thermal velocity motion caused by such waves is in- 
herent in the computed line-width variation of synthetic Si 
VIII 1445.75 A line with height in the model coronal hole, 
and that it is consistent with the observations reported by 
Banerjee et al. (1998) . Hence, we conclude that the com- 
puted spectral line broadening with height in the coronal 
hole implies the increment of non-thermal velocity as pro- 
vided by the Alfven waves. The results of our numerical 
simulations of non-linear Alfven waves show that the non- 
thermal spectral line width broadening in polar coronal holes 
can be explained by these waves. 

In our pulse-excited model, the pulse of a large ampli- 
tude is launched above the solar photosphere and this pulse 
splits into the upward and backward propagating wave- 
trains in the overlying solar atmosphere. The upward moving 
pulse train causes the instantaneous displacement of the field 
lines in perpendicular plane away and towards the line-of- 
sight. This effect generates the transversal velocity and thus 
the Doppler shift as well as line broadening. Therefore, in 
our case of pulse-excited non-linear Alfven wave model, the 
more stringent contribution of averaged transversal velocity 
V z component at a particular height causes the non-thermal 
broadening of the spectral line there. 

The theoretically estimated line-width of Si VIII 
1445.75 A synthetic line increased with height in the model 
coronal hole due to the pulse-excited, large amplitude Alfven 
waves, which is consistent with the observations of the same 
phenomenon in the polar coronal hole (Banerjee et al. 1998). 
We have also generated synthetic line profiles of Si VIII 
1445.75 A line (Fig. 4) that match well the sample obser- 
vations of Banerjee et al. (1998) showing the spectral line 
broadening. These measurements may also be the signature 
of the undamped and growing transversal pulse train that 
is the pulse-excited Alfven wave. However, the region of our 
interest is below 80 Mm in the solar atmosphere, where the 
observations show the line-width increment consistent with 
the excitation of the undamped Alfven waves. It should be 
noted that the region of decreasing line-width as reported 
by O'Shea et al. (2005) is quite high that was interpreted in 
terms of Alfven wave dissipation. We do not consider such 
regime of observations in our model, which instead is focused 
on the Alfvenic dynamics in the lower solar atmosphere of 



the coronal hole that is also the source region of the fast 
solar wind. 

It is well known that the fast solar wind starts accel- 
erating near 20 Mm in the magnetic funnels in the polar 
coronal regions (Tu et al. 2005). Therefore, the undamped 
and growing pulse-excited, large amplitude Alfven waves can 
provide the momentum to the solar wind plasma in the in- 
ner corona like the classical Alfven waves, which are typi- 
cally considered to be generated in the solar photosphere. On 
the other hand, the generation of large-amplitude Alfvenic 
pulses above the solar surface considered in this paper is 
still a puzzle. The reconnection events between the emerg- 
ing magnetic fields with the existing open magnetic field 
lines may be one of the most plausible mechanisms that can 
trigger such strong transversal velocity pulses and thus drive 
Alfven waves along the open field lines of the coronal hole. 
However, such waves and their physical properties may be 
highly dependent on the local magnetic field configuration 
of the reconnection region (Kigure et al. 2010). 

Recently, Okamoto & De Pontieu (2011) have used the 
high-resolution observations obtained by the Hinode instru- 
ments and discovered the propagating Alfven waves along 
solar spicules with a velocity amplitude of ~7 km s _1 . In 
addition, Mcintosh et al. (2011) have found evidence for the 
outward-propagating Alfvenic motions with amplitudes of 
the order of 20 km s _1 and periods of the order of 100- 
500 s throughout the quiescent atmosphere, and concluded 
that these motions carry enough energy to accelerate the 
fast solar wind and heat the quiet corona. These discoveries 
of the Alfvenic motions imply that the carried wave energy 
may cascade through various layers of the solar atmosphere. 
The main aim of this paper is to understand the role of the 
pulse-excited Alfven waves in the unresolved non-thermal 
motion in polar coronal hole and their role in the observed 
non-thermal spectral line broadening. 

In summary, the recent observations clearly show 
Alfvenic disturbances with amplitudes smaller than 50 km 
s _1 at spatial and temporal scales achieved by the simula- 
tions. This is a potentially limiting factor for the realism of 
the simulations. 

The energy carried by Alfven waves has been considered 
to be an important candidate for the heating of the coronal 
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Figure 5. The synthetic line profiles of SI VIII A1445.75 A line at three different heights y =43.5, 63.25, and 79.5 Mm in model coronal 
hole, based on the inclusion of the non-thermal contribution of Alfven wave excited by pulse amplitude A v = 25 km s -1 . The line- 
profile is deduced by considering the ionic equilibrium as reported by Mazzotta et al. (1998), the coronal hole DEM values, and coronal 
abundances as available in CHIANTI atomic data base. The simulated temperature and density are also used as input parameters. These 
synthetic line profiles does not show any significant spectral line broadening as observed by Dolla & Solomon (2008). 



hole and the acceleration of the solar wind by many authors 
(e.g., Ofman & Davila 1995, 1998; Ofman 2005, and refer- 
ences cited there). The dissipation of the energy carried by 
these waves occurs cither in the distant part of the corona 
(Parker 1991), or by some unique processes, such as phase- 
mixing, taking place in in the solar atmosphere (Heyvaerts 
& Priest 1983; Hood et al. 1997; Nakariakov et al. 1997). 
The large-amplitude non-linear Alfven waves can also be 
dissipated by coupling to longitudinal wave motions in the 
outer part of the magnetized solar corona (e.g., Boynton 
& Torkelsson 1996). This shows that the solar atmosphere 
above the polar corona may be the ideal place for the growth 
of linear and non-linear Alfven waves without much damp- 
ing, and that the waves may carry their energy to the outer 
part of the corona and can heat the solar wind ions (e.g., 
Ofman & Davila 2001). 

The main objective of the present study was to compute 
the averaged amplitude of the pulse-excited Alfven waves 
at different atmospheric heights, and explore its effects on 
the broadening of the observed spectral lines. Our compu- 
tational domain covered both the lower part of the solar at- 
mosphere as well as the inner solar corona, and we did not 
take into account any dissipative processes in our numeri- 
cal simulations because it was suggested that such processes 
may contribute to the observed narrowing of the line profiles 
(O'Shca ct al. 2005). However, this work is out of the scope 
of this paper and will be considered as a future project. 



7 CONCLUSIONS 

In this paper, we simulated numerically the behavior of im- 
pulsively excited non-linear Alfven waves in solar coronal 
holes, and studied for the first time their role in the ob- 
served broadening of spectral lines caused by these waves. 
We compared the obtained numerical results to the spec- 
tral line broadenings observed by Banerjee et al. (1998); 
Dolla & Solomon (2008). We found that the large-amplitude, 
non- linear, pulse-driven (A v —50 km s _1 ) Alfven waves are 
the most likely candidates for the non-thermal broaden- 
ing of Si VIII A1445.75 A line profiles in the polar coro- 
nal hole (Banerjee et al. 1998). Our results also show that 
Alfven waves driven by comparatively smaller velocity pulse, 



A v =25 km s _1 , only approximately contribute to the ob- 
served spectral line width of Si VIII A 1445 A at various 
heights in coronal hole as reported by Dolla & Solomon 
(2008), however, without significant evidence of broadening. 
The results of our numerical simulation and their compari- 
son to the observations of line broadening are important as 
they become an indirect evidence for the existence of larger 
amplitude, pulse-driven Alfven waves in solar coronal holes. 

Finally, we would like to suggest that more spectro- 
scopic observations should be carried out in future us- 
ing high-resolution observations (e.g., Hinodc/EIS, and also 
with upcoming Solar-C instruments) to search for a direct 
evidence of such pulse-excited non-linear Alfven waves in 
the solar atmosphere. Obviously, more detailed theoretical 
studies of a variety of pulse-driven non-linear Alfven waves 
excited by a range of pulse amplitudes are also needed in 
order to better understand the role played by these waves 
in the observed line broadening and the solar wind acceler- 
ation. 

Acknowledgments. The software used in this work 
was in part developed by the DOE-supportcd ASC/Alliancc 
Center for Astrophysical Thermonuclear Flashes at the Uni- 
versity of Chicago. This work has been supported by the 
Alexander von Humboldt Foundation (Z.E.M.) and by a 
Marie Curie International Research Staff Exchange Scheme 
Fellowship within the 7th European Community Framework 
Program (P.Ch. & K.M.) as well as by the "HPC Infras- 
tructure for Grand Challenges of Science and Engineering" 
Project, co-financed by the European Regional Develop- 
ment Fund under the Innovative Economy Operational Pro- 
gram (P.Ch. & K.M.). We also acknowledge the CHIANTI, 
which is a collaborative project involving researchers at NRL 
(USA) RAL (UK), and the Universities of Cambridge (UK), 
George Mason (USA), and Florence (Italy). A.K.S. thanks 
Shobhna Srivastava for patient encouragement, and also to 
Prof. D. Tripathi, IUCAA, Pune, India for valuable discus- 
sions on spectral line-profiles and their implications. 



REFERENCES 

Tomczyk S., MclntoshS. W., Keil S. L., et al. 2007, Science, 
317, 1192 



Nonlinear Alfven waves and the Spectral Line Broadening 9 



Antolin P., Shibata K., Kudoh T., Shiota D., & Brooks 
D. 2009, The Second Hinode Science Meeting: Beyond 
Discovery- Toward Understanding, 415, 247 

Antolin P., & Shibata K. 2010, ApJ, 712, 494 

Banerjee D., Teriaca L., Doyle J. G., & Wilhelm K. 1998, 
A&A, 339, 208 

Banerjee D., Perez-Suarez D., & Doyle J. G. 2009, A&A, 
501, L15 

Bemporad A., Abbo L., 2012, ApJ, 751, 110 
Boynton G. C, & Torkelsson U. 1996, A&A, 308, 299 
Cirtain J. W., Golub L., Lundquist L., et al. 2007, Science, 
318, 1580 

De Pontieu B., Hansteen V. H., Rouppe van der Voort L., 
van Noort M., & Carlsson M., 2007, ApJ, 655, 624 

Dere K. P., Landi E., Mason H. E., Monsignori Fossi B. C, 
& Young, P. R. 1997, A&AS, 125, 149 

Doha L., & Solomon J. 2008, A&A, 483, 271 

Dwivedi B. N., & Srivastava A. K. 2006, Sol. Phys., 237, 
143 

Dwivedi B. N., & Srivastava A. K. 2008, New A., 13, 581 
Erdelyi R., & Fedun V. 2007, Science, 318, 1572 
Hackenberg P., Marsch E., & Mann G. 2000, A&A, 360, 
1139 

Harrison R. A., Hood A. W., & Pike C. D. 2002, A&A, 
392, 319 

Hassler D. M., Rottman G. J., Shoub E. C, & Holzer T. E. 
1990, ApJ, 348, L77 

Hassler D. M., Dammasch I. E., Lemaire P., et al. 1999, 
Science, 283, 810 

Heyvaerts J., & Priest E. R. 1983, A&A, 117, 220 

Hollweg J. V. 1985, Chromospheric Diagnostics and Mod- 
elling, 235 

Hollweg J. V., & Isenberg P. A. 2007, Journal of Geophys- 
ical Research (Space Physics), 112, 8102 

Hood A. W., Ireland J., & Priest E. R. 1997, A&A, 318, 
957 

Jess D. B., Mathioudakis M., Erdelyi R., et al. 2009, Sci- 
ence, 323, 1582 

Kigure H., Takahashi K., Shibata K., Yokoyama T., & 
Nozawa, S. 2010, PASJ, 62, 993 

Kudoh T., & Shibata K. 1999, ApJ, 514, 493 

Mariska J. T., 1992, str..book 

Matsumoto T., & Shibata K. 2010, ApJ, 710, 1857 

Mazzotta P., Mazzitelli G., Colafrancesco S., Vittorio N., 
1998, A&AS, 133, 403 

McAteer R. T. J., Gallagher P. T., Williams D. R., et al. 
2003, ApJ, 587, 806 

Mcintosh S. W., de Pontieu B., Carlsson M., et al. 2011, 
Nature, 475, 477 

Moran T. G. 2003, ApJ, 598, 657 

Murawski K., & Musielak Z. E. 2010, A&A, 518, A37 
Musielak Z. E., & Moore R. L. 1995, ApJ, 452, 434 
Nakariakov V. M., Roberts B., & Murawski K. 1997, Sol. 
Phys., 175, 93 

Narain U., & Ulmschneider P. 1996, Space Sci. Rev., 75, 
453 

Ofman L., & Davila J. M. 1995, J. Geophys. Res., 100, 
23413 

Ofman L., & Davila J. M. 1998, J. Geophys. Res., 103, 
23677 

Ofman L., & Davila J. M. 2001, ApJ, 553, 935 
Ofman L. 2005, Space Sci. Rev., 120, 67 



Okamoto T. J. & De Pontieu B. 2011, ApJ, 736, L24 
O'Shea E., Banerjee D., & Doyle J. G. 2005, A&A, 436, 
L35 

Parker E. N. 1991, ApJ, 372, 719 

Pekiinhi E. R., Bozkurt Z., Afsar M., Soydugan E., & Soy- 

dugan F. 2002, MNRAS, 336, 1195 
Roberts B. 1991, Advances in Solar System Magnetohydro- 

dynamics, 105 

Roberts B. 2004, SOHO 13 Waves, Oscillations and Small- 
Scale Transients Events in the Solar Atmosphere: Joint 
View from SOHO and TRACE, 547, 1 
Srivastava A. K., Dwivedi B. N., 2007, JApA, 28, 1 
Srivastava A. K., & Dwivedi B. N. 2010, MNRAS, 405, 
2317 

Tu C.-Y., Zhou C, Marsch E., et al. 2005, Science, 308, 
519 

Verdini A., & Velli M. 2007, ApJ, 662, 669 
Verdini A., Velli M., & Buchlin E. 2009, ApJ, 700, L39 
Zaqarashvili T. V., Oliver R., & Ballester J. L. 2006, A&A, 
456, L13 



